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NUMERICAL  MODELING  STU  D I ES  over  the  past 
several  decades  have  demonstrated  progress  both  in  model  ar¬ 
chitecture  and  in  the  use  of  rapidly  advancing  computational 
resources.  Perhaps  the  most  notable  aspect  of  this  progression 
has  been  the  evolution  from  simulations  on  coarse-resolution 
horizontal  and  vertical  grids  that  outline  basins  of  simpli¬ 
fied  geometry  and  bathymetry  and  that  are  forced  by  idealized 
surface  fluxes,  to  fine- resolution  simulations  that  incorporate 
realistic  coastline  definition  and  bottom  topography  and  that 
are  forced  by  observational  data  on  relatively  short  time  scales 
(Hurlburt  and  Hogan,  2000;  Smith  et  al.,  2000;  Chassignet  and 
Garraffo,  2001;  Maltrud  and  McClean,  2005). 

The  Global  Ocean  Data  Assimilation  Experiment  (GODAE) 
is  a  coordinated  international  effort  envisioning  a  global  system 
of  observations,  communications,  modeling,  and  assimilation  that 
will  deliver  regular,  comprehensive  information  on  the  state  of  the 
oceans  in  a  way  that  will  promote  and  engender  wide  utility  and 
availability  of  this  resource  for  maximum  benefit  to  the  commu¬ 
nity.  Specific  objectives  of  GODAE  are  to  apply  state-of-the-art 
ocean  models  and  data  assimilation  methods  to  produce  short- 
range  open  ocean  forecasts,  initial  and  boundary  conditions  to 
extend  the  predictability  of  coastal  and  regional  subsystems,  and 
to  provide  initial  conditions  for  climate  forecast  models  (Inter¬ 
national  GODAE  Steering  Team,  2000).  In  the  United  States,  a 


broad  partnership  of  institutions1  is  addressing  these  objectives 
by  building  global  and  basin-scale  ocean  prediction  systems  us¬ 
ing  the  versatile  HYbrid  Coordinate  Ocean  Model  (HYCOM) 
(more  information  available  at  http://www.hycom.org). 

The  choice  of  the  vertical  coordinate  system  in  an  ocean 
model  remains  one  of  the  most  important  aspects  of  its  design. 
In  practice,  the  representation  and  parameterization  of  the  pro¬ 
cesses  not  resolved  by  the  model  grid  are  often  directly  linked 
to  the  vertical  coordinate  choice  (Griffies  et  al.,  2000).  Oceanic 
general  circulation  models  traditionally  represent  the  vertical 
in  a  series  of  discrete  intervals  in  either  a  depth,  density,  or  ter¬ 
rain-following  unit.  Recent  model  comparison  exercises  per¬ 
formed  in  Europe  (Willebrand  et  al.,  2001)  and  in  the  United 
States  (Chassignet  et  al.,  2000)  have,  however,  shown  that  the 
use  of  only  one  vertical  coordinate  representation  cannot  be 
optimal  everywhere  in  the  ocean.  These  and  earlier  compari¬ 
son  studies  (Chassignet  et  al.,  1996;  Marsh  et  al.,  1996,  Roberts 
et  al.,  1996)  show  that  all  the  models  considered  were  able  to 
simulate  large-scale  characteristics  of  the  oceanic  circulation 
reasonably  well,  but  the  interior  water-mass  distribution  and 
associated  thermohaline  circulation  are  strongly  influenced 
by  localized  processes  that  are  not  represented  equally  by  each 
model’s  vertical  discretization. 


1  University  of  Miami,  Naval  Research  Laboratory  (NRL),  National  Aeronautics  and  Space  Administration/Goddard  Institute  for  Space  Studies  (NASA/GISS),  National  Oceanic  and  Atmospheric 
Administration/National  Centers  for  Environmental  Prediction  (NOAA/NCEP),  National  Oceanic  and  Atmospheric  Administration/Atlantic  Oceanographic  Meteorological  Laboratory 
(NOAA/AOML),  National  Oceanic  and  Atmospheric  Administration/Pacific  Marine  Environmental  Laboratory  (NOAA/PMEL),  Planning  Systems,  Inc.  (PSI),  Fleet  Numerical  Meteorology  and 
Oceanography  Center  (FNMOC),  Naval  Oceanographic  Office  (NAVOCEANO),  Service  Hydrographique  et  Oc£anographique  de  la  Marine  (SHOM),  Laboratoire  des  Ecoulements  Geophy¬ 
siques  et  Industrie^  (LEGI),  Open  Source  Project  for  a  Network  Data  Access  Protocol  (OPeNDAP),  University  of  North  Carolina,  Rutgers,  University  of  South  Florida,  Fugro  GEOS,  Roffer's  Ocean 
Fishing  and  Forecasting  Service  (ROFFS),  Orbimage,  Shell,  ExxonMobil. 
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Because  none  of  the  three  main  verti¬ 
cal  coordinates  (depth,  density,  and  ter¬ 
rain-following)  (see  Box  1  for  details) 
provide  universal  optimality,  it  is  natural 
to  envision  a  hybrid  approach  that  com¬ 
bines  the  best  features  of  each  vertical 
coordinate.  Isopycnic  (density-track¬ 
ing)  layers  work  best  for  modeling  the 
deep  stratified  ocean,  levels  at  constant 
fixed  depth  or  pressure  are  best  to  use  to 
provide  high  vertical  resolution  near  the 
surface  within  the  mixed  layer,  and  ter¬ 
rain-following  levels  are  often  the  best 
choice  for  modeling  shallow  coastal  re¬ 
gions.  In  HYCOM,  the  optimal  vertical 
coordinate  distribution  of  the  three  ver¬ 
tical  coordinate  types  is  chosen  at  every 
time  step.  The  hybrid  vertical  coordinate 
generator  makes  a  dynamically  smooth 
transition  among  the  coordinate  types 
using  the  continuity  equation. 


HYBRID  VERTICAL 
COORDINATES 

Hybrid  vertical  coordinates  can  mean 
different  things  to  different  people:  they 
can  be  a  linear  combination  of  two  or 
more  conventional  coordinates  (Song 
and  Haidvogel,  1994;  Ezer  and  Mellor, 
2004;  Barron  et  al.,  2006)  or  they  can  be 
truly  generalized  (i.e.,  aiming  to  mimic 
different  types  of  coordinates  in  different 
regions  of  a  model  domain)  (Bleck,  2002; 
Burchard  and  Beckers,  2004;  Adcroft  and 
Hallberg,  2006;  Song  and  Hou,  2006). 

The  generalized  vertical  coordinates  in 
HYCOM  deviate  from  isopycnals  (con¬ 
stant  density  surfaces)  wherever  the  latter 
may  fold,  outcrop,  or  generally  provide 
inadequate  vertical  resolution  in  portions 
of  the  model  domain.  HYCOM  is  at  its 
core  a  Lagrangian  layer  model,  except 
for  the  remapping  of  the  vertical  coor¬ 


dinate  by  the  hybrid  coordinate  genera¬ 
tor  after  all  equations  are  solved  (Bleck, 
2002;  Chassignet  et  al.,  2003;  Halliwell, 
2004)  and  for  the  fact  that  there  is  a  non¬ 
zero  horizontal  density  gradient  within 
all  layers.  HYCOM  is  thus  classified  as 
a  Lagrangian  Vertical  Direction  (LVD) 
model  in  which  the  continuity  (thickness 
tendency)  equation  is  solved  forward  in 
time  throughout  the  domain,  while  an 
Arbitrary  Lagrangian- Eulerian  (ALE) 
technique  is  used  to  re-map  the  vertical 
coordinate  and  maintain  different  co¬ 
ordinate  types  within  the  domain.  This 
differs  from  Eulerian  Vertical  Direction 
(EVD)  models  with  fixed  vertical  coor¬ 
dinates  that  use  the  continuity  equation 
to  diagnose  vertical  velocity  (Adcroft  and 
Hallberg,  2006).  The  ability  to  adjust  the 
vertical  spacing  of  the  coordinate  sur¬ 
faces  in  HYCOM  simplifies  the  numeri- 


BOX  1:  OCEAN  REGIMES  AND  VERTICAL  COORDINATES  [ 


Schematic  of  an  ocean  basin  illustrating 
the  three  regimes  of  the  ocean  germane 
to  the  considerations  of  an  appropriate 
vertical  coordinate.  The  surface  mixed 
layer  is  naturally  represented  using  fixed- 
depth  z  (or  pressure  p)  coordinates,  the 
interior  is  naturally  represented  using 
isopycnic  p  (density  tracking)  coordi¬ 
nates;  and  the  bottom  boundary  is  natu¬ 
rally  represented  using  terrain-following 
a  coordinates  (after  Griffies  et  al.,  2000). 
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cal  implementation  of  several  physical 
processes  (e.g.,  mixed-layer  detrainment, 
convective  adjustment,  sea- ice  modeling) 
without  harming  the  model  of  the  basic 
and  numerically  efficient  resolution  of 
the  vertical  that  is  characteristic  of  iso- 
pycnic  models  throughout  most  of  the 
oceans  volume  (Bleck  and  Chassignet, 
1994;  Chassignet  et  al.,  1996). 

The  default  configuration  of  HYCOM 
is  isopycnic  in  the  open  stratified  ocean, 
but  makes  a  dynamically  smooth  tran¬ 
sition  to  terrain- following  coordinates 
in  shallow  coastal  regions  and  to  fixed 
pressure-level  coordinates  in  the  surface 
mixed  layer  and/or  unstratified  seas  (Fig¬ 
ure  1).  In  doing  so,  the  model  takes  ad¬ 
vantage  of  the  different  coordinate  types 
in  optimally  simulating  coastal  and  open- 
ocean  circulation  features.  A  user-chosen 
option  allows  specification  of  the  vertical 
coordinate  separation  that  controls  the 
transition  among  the  three  coordinate 
systems.  The  assignment  of  additional 
coordinate  surfaces  to  the  oceanic  mixed 
layer  also  allows  the  straightforward  im¬ 
plementation  of  multiple  vertical  mixing 
turbulence  closure  schemes  (Halliwell, 
2004).  The  choice  of  the  vertical  mixing 
parameterization  is  also  of  importance 
in  areas  of  strong  entrainment,  such  as 
overflows  (Xu  et  al.,  submitted). 

Figure  1  illustrates  the  transition 
among  pressure,  terrain -following,  and 
isopycnic  coordinates  in  1/25°  West 
Florida  Shelf  simulations  nested  within 
a  1/12°  North  Atlantic  configuration, 
and  it  demonstrates  the  flexibility  with 
which  vertical  coordinates  can  be  chosen 
and  the  capability  of  adding  additional 
vertical  resolution.  The  original  vertical 
discretization  used  in  the  1/12°  North 
Atlantic  configuration  is  compared  to 


two  others  with  six  layers  added  at  the 
top — one  with  pressure-level  coordinates 
and  the  other  with  terrain-following 
coordinates  over  the  shelf — because  in 
many  coastal  applications  it  is  desirable 
to  provide  higher  resolution  from  sur¬ 
face  to  bottom  to  adequately  resolve  the 
vertical  structure  of  water  properties  and 
of  the  bottom  boundary  layer  in  shallow 
water.  Halliwell  et  al.  (in  preparation) 
document  the  advantages  and  disadvan¬ 
tages  of  the  choices  shown  in  Figure  1 . 

COMPUTATIONAL  ASPECTS 

Before  one  can  numerically  solve  (i.e., 
discretize)  the  models  equations,  one 
must  decide  on  a  global  projection  and 
on  how  to  treat  the  singularity  associ¬ 
ated  with  the  North  Pole  (the  South 
Pole,  being  over  land,  is  not  an  issue). 
There  are  several  projections  that  al¬ 


low  the  Arctic  to  be  included  in  a  global 
ocean  model.  The  HYCOM  global  con¬ 
figuration  uses  an  Arctic  dipole  patch 
matched  to  a  standard  Mercator  grid 
at  47°N  (Figure  2).  Because  our  tar¬ 
get  horizontal  resolution  is  1/1 2°  at  the 
equator  (i.e.,  9  km),  the  location  of  the 
dipoles  at  47°N  gives  us  a  good  resolu¬ 
tion  at  mid  latitude  (i.e.,  7  km)  as  well 
as  in  the  Arctic  Ocean  (i.e.,  3.5  km  at  the 
North  Pole)  where  the  Rossby  radius  of 
deformation  is  smaller.  An  advantage 
of  this  pole-shifting  projection,  as  op¬ 
posed  to  most  others,  is  that  all  the  grid 
points  below  47°N  remain  on  regular 
grid  spacing  (i.e.,  x-y  grid  ratio  of  1). 
The  corresponding  numerical  array  size 
is  4500  by  3298  with  32  hybrid  layers  in 
the  vertical.  The  complete  system  will 
include  the  Los  Alamos  Sea  Ice  Model 
(CICE)  (Hunke  and  Lipscomb,  2004)  on 
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Figure  1.  Cross  sections  of  layer  density  and  model  interfaces  across  the  West  Florida  Shelf 
in  a  1/25°  West  Florida  Shelf  subdomain  covering  the  Gulf  of  Mexico  east  of  87°W  and 
north  of  23°N  and  embedded  in  a  1/12°  Atlantic  basin  HYCOM  simulation  (Halliwell  et  al., 
in  preparation).  This  figure  illustrates  the  transition  among  pressure,  terrain-following,  and 
isopycnic  coordinates  in  1/25°  West  Florida  Shelf  simulations  nested  within  a  1/12°  North 
Atlantic  configuration,  and  it  demonstrates  the  flexibility  with  which  vertical  coordinates 
can  be  chosen  and  the  capability  of  adding  additional  vertical  resolution. 


the  same  grid.  The  ocean  and  ice  models 
will  run  simultaneously,  but  on  sepa¬ 
rate  sets  of  processors  (a  much  smaller 
number  for  CICE),  communicating  via 
an  Earth  System  Modeling  Framework 
(ESMF) -based  coupler  (Hill  et  al.,  2004). 

HYCOM’s  basic  parallelization  strat¬ 
egy  is  two-dimensional  domain  decom¬ 
position  (i.e.,  the  entire  model  domain 
is  divided  up  into  smaller  sub-domains, 
or  tiles,  and  each  processor  “owns”  one 
tile).  A  halo  region  is  added  around  each 
tile  to  allow  communication  operations 
(e.g.,  updating  the  halo)  to  be  com¬ 
pletely  separated  from  computational 
kernels,  greatly  increasing  the  maintain¬ 
ability  and  expandability  of  the  code 
base.  Rather  than  the  conventional  one- 
or  two-grid  points- wide  halo,  HYCOM 
has  a  six-grid-point-wide  halo,  which  is 
“consumed”  over  several  operations  to 
reduce  halo  communication  overhead. 
For  global  and  basin-scale  applications, 
it  is  important  to  avoid  calculations  over 
land.  HYCOM  fully  “shrink  wraps”  cal¬ 
culations  on  each  tile  and  discards  tiles 
that  are  completely  over  land  (Bleck  et 
al.,  1995).  HYCOM  goes  farther  than 
most  structured  grid  ocean  models  in 
land  avoidance  by  allowing  more  than 
one  neighboring  tile  to  the  north  and 
south.  Figure  2  shows  the  current  tiling 
for  the  1/12°  global  domain  with  equal 
sized  tiles  that  (a)  allows  rows  to  be  offset 
from  each  other  if  this  gives  fewer  tiles 
over  the  ocean  and  (b)  allows  two  tiles  to 
be  merged  into  one  larger  tile  if  less  than 
50%  of  their  combined  area  is  ocean.  The 
Mediterranean  region  of  Figure  2  illus¬ 
trates  both  these  optimizations. 

In  the  example  presented  in  Figure  3, 
the  global  model  is  initialized  from  an 
oceanic  climatology  of  temperature 
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and  salinity.  The  model  is  forced  by 
prescribed  climatological  atmospheric 
wind,  thermal,  and  precipitation  fields 
while  evaporation  is  computed  using 
the  modeled  sea  surface  temperature 
(SST).  There  is  also  a  weak  relaxation 
to  climatological  surface  salinity.  These 
simulations  use  a  simple  thermodynamic 
sea-ice  model  (CICE  is  running  stand¬ 
alone  on  the  global  Arctic  patch  grid,  but 
we  are  awaiting  ESMF-based  coupling 
between  HYCOM  and  CICE  before  run¬ 
ning  a  coupled  case).  Even  without  ice 
dynamics,  the  seasonal  cycle  of  ice  cover¬ 
age  is  good  overall  (not  illustrated).  One 
reason  for  the  good  agreement  is  that  the 
atmospheric  forcing  is  based  on  an  ac¬ 
curate  ice  extent,  which  provides  a  strong 


tendency  for  the  ocean/sea- ice  system 
to  form  ice  appropriately.  Figure  3 
compares  the  sea  surface  height  (SSH) 
variability  from  the  climatologically 
forced  1/12°  global  HYCOM  to  the  Oct. 
1992-Nov.  1998  SSH  variability  based  on 
Topex-Poseidony  ERS-1 ,  and  ERS-2  altim¬ 
eter  data  (derived  by  Collecte  Localisa¬ 
tion  Satellite  (CLS),  France).  Overall,  the 
modeled  regions  of  high  variability  are 
in  reasonably  good  agreement  with  the 
observations,  especially  in  the  Antarctic 
Circumpolar  Current  region.  The  equa¬ 
torial  Pacific  is  an  exception  because  the 
altimeter  data  include  interannual  vari¬ 
ability  not  present  in  the  model  (e.g., 
the  large  variability  associated  with  the 
1997-1998  El  Nino). 


OCEAN  PREDICTION 

Although  HYCOM  is  a  relatively  sophis¬ 
ticated  model  that  includes  a  large  suite 
of  physical  processes  and  incorporates 
numerical  techniques  that  are  optimal 
for  dynamically  different  regions  of  the 
ocean,  data  assimilation  is  still  essential 
for  ocean  prediction  because  (a)  many 
ocean  phenomena  are  due  to  nonlinear 
processes  (i.e.,  flow  instabilities)  and  thus 
are  not  a  deterministic  response  to  at¬ 
mospheric  forcing,  (b)  errors  exist  in  the 
atmospheric  forcing,  and  (c)  ocean  mod¬ 
els  are  imperfect,  including  limitations  in 
numerical  algorithms  and  in  resolution. 

Substantial  information  about  the 
ocean  surface's  space-time  variability 
is  obtained  remotely  from  instruments 


Figure  2.  Current  tiling  for 
the  1/12°  global  domain 
with  equal-sized  tiles 
that  (a)  allows  rows  to  be 
offset  from  each  other  if 
this  gives  fewer  tiles  over 
the  ocean  and  (b)  allows 
two  tiles  to  be  merged 
into  one  larger  tile  if  less 
than  50  percent  of  their 
combined  area  is  ocean. 
The  Mediterranean  region 
especially  illustrates  both 
these  optimizations;  out 
of  the  original  1152  (36  by 
32)  approximately  equal¬ 
sized  tiles,  371  tiles  are 
entirely  land  and  are  dis¬ 
carded,  leaving  781.  The 
Arctic  “wraps"  between 
the  left  and  right  halves  of 
the  top  edge. 
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Figure  3.  Comparison  be¬ 
tween  the  (observed)  Oct. 
1992-May  2005  sea  surface 
height  (SSH)  variability 
based  on  Topex-Poseidon, 
ERS-1,  and  ERS-2  altimeter 
data  (top)  (derived  by  Col- 
lecte  Localisation  Satellite 
[CLS],  France)  and  three 
years  of  (modeled)  sea  sur¬ 
face  height  (SSH)  variability 
from  the  climatologically 
forced  1/1 2°  global  HYCOM 
(bottom).  Overall,  the 
modeled  regions  of  high 
variability  are  in  good  agree¬ 
ment  with  the  observations, 
especially  in  the  Antarctic 
Circumpolar  Current  region. 
The  equatorial  Pacific  is  an 
exception  because  the  altim¬ 
eter  data  include  interannual 
variability  not  present  in  the 
model,  for  example,  the  large 
variability  associated  with 
the  1997- 1998  El  Nino. 


aboard  satellites,  but  these  observations 
are  insufficient  for  specifying  the  sub¬ 
surface  variability.  Vertical  profiles  from 
expendable  bathythermographs  (XBT), 
conductivity-temperature-depth  (CTD) 
profilers,  and  profiling  floats  (e.g.,  Argo, 
which  measures  temperature  and  salinity 
in  the  upper  2000  m  of  the  ocean)  pro¬ 
vide  another  substantial  source  of  data. 
Even  together,  these  data  sets  are  insuffi¬ 
cient  to  determine  the  state  of  the  ocean 
completely,  so  it  is  necessary  to  exploit 
prior  knowledge  in  the  form  of  statis¬ 


tics  determined  from  past  observations 
as  well  as  our  present  understanding  of 
ocean  dynamics.  By  combining  all  of 
these  observations  through  data  assimi¬ 
lation  into  an  ocean  model,  it  is  possible 
to  produce  a  dynamically  consistent  de¬ 
piction  of  the  ocean.  It  is,  however,  ex¬ 
tremely  important  that  the  freely  evolv¬ 
ing  ocean  model  (i.e.,  non-data-assimila- 
tive  model)  has  skill  in  hindcasting  and 
in  predicting  ocean  features  of  interest. 
For  a  detailed  overview  of  the  HYCOM 
data  assimilative  system,  the  reader  is 


referred  to  Chassignet  et  al.  (in  press). 

The  present  Navy  near- real -time 
1/12°  North  Atlantic  HYCOM  ocean 
forecasting  system  is  the  first  step  toward 
the  fully  global  1/12°  HYCOM  predic¬ 
tion  system.  The  North  Atlantic  system 
assimilates  daily,  real-time  satellite  al¬ 
timeter  data  (Geosat-Follow-On  [GFO], 
Environmental  Satellite  [ENVISAT]  i, 
and  Jason- 1).  These  data  are  provided 
via  the  Altimeter  Data  Fusion  Center 
(ADFC)  at  the  Navel  Oceanographic 
Office  (NAVOCEANO)  to  generate  the 
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two-dimensional  Modular  Ocean  Data 
Assimilation  System  (MODAS)  SSH 
analysis  (Fox  et  al.,  2002).  The  MODAS 
analysis  is  an  optimal  interpolation 
technique  that  uses  complex  covariance 
functions,  including  spatially  varying 
length  and  time  scales  as  well  as  propa¬ 
gation  terms  derived  from  many  years  of 
altimetry  (Jacobs  et  al.,  2001).  The  mod¬ 
el  SST  is  relaxed  to  the  daily  MODAS 
SST  analysis  that  uses  daily  Multi-Chan¬ 
nel  Sea  Surface  Temperature  (MCSST) 
data  derived  from  the  5-channel  Ad¬ 
vanced  Very  High  Resolution  Radiome¬ 
ters  (AVHRR) — globally  at  8.8  km  reso¬ 
lution  and  at  2  km  in  selected  regions. 

To  properly  assimilate  the  SSH 
anomalies  determined  from  satellite 
altimeter  data,  the  oceanic  mean  SSH 
over  the  altimeter  observation  period 
must  be  determined.  In  this  mean,  it  is 
essential  that  the  mean  current  systems 
and  associated  SSH  fronts  be  accurately 
represented  (position,  amplitude,  and 
sharpness).  Unfortunately,  Earths 
geoid  is  not  presently  known  with 
sufficient  accuracy  for  this  purpose, 
and  coarse  hydrographic  climatologies 
(~1°  horizontal  resolution)  cannot 
provide  the  spatial  resolution  necessary. 
HYCOM,  therefore,  uses  a  mean  SSH 
from  a  previous  fully  eddy- resolving 
ocean  model  simulation  that  was  found 
to  have  fronts  in  the  correct  position 
and  is  consistent  with  hydrographic 
climatologies  (Chassignet  and  Garraffo, 
2001).  Several  satellite  missions  are 
either  underway  or  planned  to  determine 
a  more  accurate  geoid,  but  until  the 
accuracy  reaches  a  few  centimeters  on 
horizontal  scales  of  about  30  km,  the 
present  approach  will  be  necessary. 

The  North  Atlantic  system  runs  weekly 


every  Wednesday  and  produces  a  10-day 
hindcast  and  a  14-day  forecast.  The  Navy 
Operational  Global  Atmospheric  Predic¬ 
tion  System  (NOGAPS)  (Rosmond  et  al., 
2002)  provides  the  atmospheric  forcing, 
but  in  the  14-day  forecasts,  the  forcing 
linearly  reverts  toward  climatology  af¬ 
ter  five  days.  During  the  forecast  period, 
the  SST  is  relaxed  toward  climatologi- 
cally  corrected  persistence  of  the  nowcast 
SST  with  a  relaxation  time  scale  of  one- 
fourth  the  forecast  length  (i.e.,  one  day 
for  a  four-day  forecast).  The  impact  of 
these  choices  is  discussed  by  Smedstad  et 
al.  (2003)  and  Shriver  et  al.  (in  press).  For 
an  evaluation  of  the  North  Atlantic  sys¬ 
tem,  the  reader  is  referred  to  Chassignet 
et  al.  (2005,  in  press). 

The  near  real-time  North  Atlantic  ba¬ 
sin  model  outputs  are  made  available  to 
the  ocean  science  community  within  24 
hours  via  the  HYCOM  Consortium  data 
server  (more  information  available  at 
http://www.hycom.org/dataserver)  using 
a  familiar  set  of  tools  such  as  OPeNDAP, 
Live  Access  Server  (LAS),  and  file  trans¬ 
fer  protocol  (FTP).  These  tools  have 
been  modified  to  perform  with  hybrid 
vertical  coordinates  to  provide  HYCOM 
subsets  to  coastal  or  regional  nowcast/ 
forecast  partners  as  initial  and  boundary 
conditions.  The  LAS  has  been  imple¬ 
mented  with  an  intuitive  user  interface 
to  enhance  the  usability  of  ocean  pre¬ 
diction  system  outputs  and  to  perform 
diagnostics. 

BOUNDARY  CONDITIONS  FOR 
REGIONAL  MODELS 

The  chosen  horizontal  and  vertical  reso¬ 
lution  for  the  above  HYCOM  predic¬ 
tion  system  only  marginally  resolves  the 
coastal  ocean  (7  km  at  mid  latitudes, 


with  up  to  15  terrain-following  a  coordi¬ 
nates  over  the  shelf),  but  provides  an  ex¬ 
cellent  starting  point  for  even  higher- res¬ 
olution  coastal  ocean  prediction  systems. 
The  model  resolution  should  increase 
to  1/25°  (3-4  km  at  mid-latitudes)  by 
the  end  of  the  decade.  An  important  at¬ 
tribute  of  the  data  assimilative  HYCOM 
simulations  is,  therefore,  its  capability  to 
provide  boundary  conditions  to  regional 
and  coastal  models. 

To  increase  the  predictability  of 
coastal  regimes,  several  partners  within 
the  HYCOM  consortium  are  develop¬ 
ing  and  evaluating  boundary  conditions 
for  coastal  prediction  models  based  on 
the  HYCOM  data  assimilative  system 
outputs.  The  inner  nested  models  may 
or  may  not  be  HYCOM,  so  the  coupling 
of  the  global  and  coastal  models  must  be 
able  to  handle  dissimilar  vertical  grids. 
Coupling  HYCOM  to  HYCOM  is  now 
routine  via  one-way  nesting  (Zamu- 
dio  et  al.,  in  preparation).  Outer  model 
fields  are  periodically  interpolated  to 
the  horizontal  mesh  of  the  nested  model 
(specified  by  the  user,  but  typically  daily) 
and  stored  in  an  archive  file.  The  num¬ 
ber  of  coordinates  can  be  increased  to 
augment  the  vertical  resolution  of  the 
nested  model  and  to  ensure  that  there  is 
sufficient  vertical  resolution  to  resolve 
the  bottom  boundary  layer.  The  nested 
model  is  initialized  from  the  first  archive 
file;  the  entire  set  of  archives  provides 
boundary  conditions  during  the  nested 
run,  ensuring  consistency  between  initial 
and  boundary  conditions.  Coupling  HY¬ 
COM  to  other  finite  difference  models, 
such  as  the  Navy  Coastal  Ocean  Model 
(NCOM)  or  the  Regional  Ocean  Model 
System  (ROMS),  has  already  been  dem¬ 
onstrated,  and  coupling  of  HYCOM  to 
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unstructured  grid/finite  element  models 
is  in  progress. 

We  now  describe  the  use  of  near- real¬ 
time  HYCOM  nowcasts  and  forecasts  as 
boundary  and  initial-condition  provid¬ 
ers  to  a  nested  coastal  simulation  in  the 
South  Atlantic  Bight  (SAB)  region  of 
the  eastern  U.S.  coast.  The  1/12°  North 
Atlantic  HYCOM  does  not  necessarily 
include  all  forcing  and  physics  relevant 
at  the  coastal  region  scale  (e.g.,  lack  of 
tidal  forcing  in  the  present  simulation). 
The  nesting  of  higher-resolution  models 
within  the  basin-scale  HYCOM  there¬ 
fore  allows  limited-area  regional  forc¬ 
ings  (terrestrial  buoyancy  inputs,  tides), 
physics  (wetting  and  drying),  and  coastal 
geometry  (tidal  inlets,  estuaries)  to  add 
value  to  the  larger-scale  HYCOM  ocean- 
state  estimates. 

The  quasi-operational  regional-scale 
modeling  system  developed  at  the  Uni¬ 


versity  of  North  Carolina  (UNC)-SAB 
(National  Ocean  Partnership  Program 
[NOPP] -funded  South  Atlantic  Bight 
Limited  Area  Model  [SABLAM],  South- 
East  U.S.  Atlantic  Coastal  Ocean  Observ¬ 
ing  System  [SEACOOS])  (Blanton,  2003) 
uses  the  finite  element  coastal  ocean 
model  QUODDY  (Lynch  et  al.,  1996). 
Terrestrial  buoyancy  inputs  to  the  conti¬ 
nental  shelf,  strong  tides,  and  a  vigorous 
western  boundary  current  contribute  to 
the  complexity  of  this  region.  To  simu¬ 
late  the  density-dependent  dynamics, 
the  UNC-SAB  modeling  system  is  nested 
within  the  HYCOM  GODAE  near-real¬ 
time  system.  Boundary  and  initialization 
data  for  the  SAB  regional-scale  model 
are  obtained  from  the  HYCOM  GODAE 
Live  Access  Server  and  are  mapped  to  the 
finite  element  regional  model  domain. 
For  each  forecast,  the  system  spins  up  the 
regional  tides  for  five  model  days,  with 


the  density  field  held  fixed.  The  atmo¬ 
spheric  (buoyancy  and  momentum)  and 
river  fluxes  are  turned  on.  The  timing  is 
synchronized  such  that  the  fluxes  are  ac¬ 
tive  for  one  day  at  the  time  the  HYCOM 
fields  are  valid. 

An  example  of  this  nested  system  is 
shown  in  Figure  4.  The  HYCOM  now- 
cast  for  November  28,  2005  is  used  to 
initialize  the  regional  model  domain, 
onto  which  the  tides  and  high-resolu¬ 
tion  atmospheric  fluxes  are  applied.  The 
resulting  solution  is  used  to  initialize  and 
drive  the  limited-area,  estuary- resolving 
finite  element  implementation.  The  ef¬ 
fects  of  both  the  Gulf  Stream,  as  provid¬ 
ed  by  the  HYCOM  initial  and  boundary 
conditions,  and  the  local  tides  are  seen  in 
the  limited-area  model  (Figure  4,  right). 
Strong  along-shelf  and  poleward  flow 
is  seen  at  the  shelfbreak.  The  poleward, 
offshelf- directed  flow  on  the  shelf  is  due 


Figure  4.  (Left)  1/12°  North  Atlantic  HYCOM-GODAE  sea  surface  temperature  (SST)  nowcast  for  November  28,  2005.  The  three-dimensional  solution  for 
this  date  is  used  to  initialize  the  University  of  North  Carolina-South  Atlantic  Bight  (UNC-SAB)  finite  element  modeling  system  (middle,  every  5th  model 
vector).  The  surface  temperature  and  velocity  are  shown,  after  addition  of  the  regional  tides  and  atmospheric  fluxes.  The  limited-area  finite  element  imple¬ 
mentation  (right,  every  3rd  model  vector)  includes  the  estuary  and  tidal  inlets  along  the  Georgia/South  Carolina  coast  and  extends  to  the  shelf-break.  This 
snapshot  also  shows  surface  temperature  and  velocity  and  is  for  three  days  after  initialization.  The  color  scales  between  the  HYCOM  and  nested  plots  are 
not  the  same. 
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to  the  tides. 

Evaluation  of  the  near  real-time  HY- 
COM  outputs  relative  to  available  ob¬ 
servations  in  the  SAB  consists  of  com¬ 
parisons  to  National  Ocean  Service  water 
levels  along  the  coast  and  mid-shelf  tem¬ 
perature  from  an  established  continental 
shelf  observational  network  (SABSOON) 
(Seim,  2000).  Variability  in  observed 
coastal  water  levels  is  due  to  tides,  wind- 
stress-driven  and  other  lower- frequency 
fluctuations  (deep-ocean  contributions 
to  shelf- wide  sea  level).  In  the  SAB,  tides 
account  for  at  least  90%  of  the  total  water 
level  variability.  Figure  5a  shows  subtidal 
coastal  water  levels  for  two  stations  in 
the  SAB.  Lower-frequency,  seasonal -scale 
variations  are  well  captured  by  HYCOM. 
Weather- band  fluctuations  are  also  rea¬ 
sonably  well  represented.  Observations  of 
in  situ  water  fields  (salinity,  temperature) 
are  comparatively  less  available  in  the 


SAB.  The  SABSOON  observational  net¬ 
work,  situated  on  the  Georgia  continental 
shelf,  has  been  making  routine  and  real¬ 
time  observations  of  water  properties  for 
three  years.  Figure  5b  shows  observed 
SABSOON  R2  near-surface  temperature 
and  HYCOM  mixed-layer  temperature. 
The  HYCOM  mixed  layer  covers  most 
of  the  water  column  vertical  grid  in  this 
region.  The  root  mean  square  (RMS)  er¬ 
ror  between  the  signals  (where  both  are 
available)  is  1.4°  C.  Strong,  prolonged 
cooling  of  the  SAB  continental  shelf  dur¬ 
ing  summer  2003  is  seen  in  the  HYCOM 
temperature  (Figure  5b),  and  is  attrib¬ 
utable  to  a  variety  of  coincident  envi¬ 
ronmental  conditions  (for  a  review  see 
Aretxabaleta  et  al.  [submitted]).  The  HY¬ 
COM  best-prior-estimates  from  summer 
2003  are  being  used,  with  the  regional 
modeling  system,  to  examine  the  physical 
nature  of  this  extreme  cooling  event. 


OUTLOOK 

The  long-term  goal  of  the  HYCOM 
consortium  is  an  eddy- resolving,  fully 
global  ocean  prediction  system  with 
data  assimilation  to  be  transitioned  to 
the  U.S.  Naval  Oceanographic  Office  at 
1/12°  equatorial  (~7  km  mid-latitude) 
resolution  in  2007  and  1/25°  resolution 
by  201 1.  Development  of  the  global  sys¬ 
tem  is  underway  and  includes  the  ocean 
model,  ice  model,  tides,  and  data  assimi¬ 
lation.  Data  assimilation  is  traditionally 
formulated  as  a  least-squares  estimation 
problem.  In  spite  of  a  fairly  simple  theo¬ 
retical  framework,  application  to  non¬ 
linear  numerical  models  of  the  ocean 
circulation  is  far  from  trivial  (Brasseur, 
2006).  The  difficulty  is  in  finding  algo¬ 
rithms  that  provide  an  acceptable  solu¬ 
tion  in  terms  of  computer  resources.  The 
size  of  the  problem  makes  it  indeed  very 
difficult  to  use  sophisticated  assimilation 
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Figure  5.  (a)  Water  level  comparison  for  two  stations  in  the  SAB.  Daily  HYCOM  best-estimate  water  levels  (blue)  and  observed  NOS  water  levels  (red) 
are  shown  for  Fort  Pulaski,  Georgia  and  Virginia  Key,  Florida.  The  two  stations  are  arbitrarily  offset  for  clarity,  (b)  Observed,  near-surface  temperature 
(blue)  and  HYCOM  mixed-layer  temperature  (red)  for  the  SABSOON  mid-shelf  station  R2.  The  root  mean  square  (RMS)  error  at  this  location  is  1.4°C. 
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techniques  because  some  of  these  meth¬ 
ods  can  increase  the  cost  of  running  the 
model  by  a  factor  of  100.  The  strategy  we 
adopted  is  to  start  with  a  low-cost  and 
simple  data  assimilation  approach  (e.g., 
Cooper  and  Haines,  1996),  and  then 
gradually  increase  the  complexity.  Sev¬ 
eral  sophisticated  data  assimilation  tech¬ 
niques  are  already  in  place  to  work  with 
HYCOM  and  are  being  evaluated.  These 
techniques  are,  in  increasing  level  of  so¬ 
phistication,  the  Naval  Research  Labo¬ 
ratory  (NRL)’s  Coupled  Ocean  Data 
Assimilation  (NCODA),  the  Singular 
Evolutive  Extended  Kalman  (SEEK)  fil¬ 
ter,  the  Reduced  Order  Information  Fil¬ 
ter  (ROIF),  the  Reduced  Order  Adaptive 
Filter  (ROAF)  (including  adjoint),  the 
Ensemble  Kalman  Filter  (EnKF),  and  the 
4D-VAR  Representer  method.  Although 
these  techniques  work  with  HYCOM,  it 
does  not  mean  that  they  will  be  used  op¬ 
erationally:  the  NCODA  and  SEEK  tech¬ 
niques  are  presently  being  considered  as 
the  next  generation  data  assimilation  to 
be  used  in  the  near- real-time  system.  The 
remaining  techniques,  because  of  their 
cost,  are  being  evaluated  mostly  within 
specific  limited  areas  of  high  interest  or 
coastal  HYCOM  configurations. 

Another  HYCOM  North  Atlantic  con¬ 
figuration  forms  the  backbone  of  the 
National  Oceanic  and  Atmospheric  Ad¬ 
ministration  (NOAA)/National  Centers 
for  Environmental  Prediction  (NCEP)/ 
North  Atlantic  Ocean  Forecast  System 
(NAOFS)  (more  information  available  at 
http://polar.ncep.noaa.gov/ofs/),  which 
provides  daily  nowcasts  and  five-day 
forecasts.  It  primarily  differs  from  the 
North  Atlantic  Navy  system  described 
earlier  in  three  ways:  (a)  the  choice  of 
horizontal  grid,  (b)  the  use  of  NCEP- 


based  wind  and  thermal  forcing,  and  (c) 
the  choice  of  data  assimilation  technique. 
The  NOAA/NCEP  group  is  using  a  con¬ 
figuration  that,  for  the  same  number 
of  grid  points  as  in  the  regular  Merca¬ 
tor  projection  used  in  the  Navy  system, 
has  finer  resolution  in  the  western  and 
northern  portions  of  the  basin  and  on 
shelves  (3-7  km),  in  order  to  provide 
higher  resolution  along  the  U.S.  coast 
rather  than  toward  the  east  and  southeast 
(7-13  km).  The  model  domain  configu¬ 
ration  is  from  20°S  to  76°N,  including 
marginal  seas,  except  for  the  Mediterra¬ 
nean  and  Baltic  Seas.  Atmospheric  mo¬ 
mentum,  heat,  and  water  fluxes  are  de¬ 
rived  from  the  three  hourly  NCEP-based 
fields.  Tidal  forcing  and  river  outflows 
are  prescribed.  The  observations  used  in 
the  assimilation  at  present  are  remotely 
sensed  and  in  situ  SST,  and  work  is  un¬ 
derway  to  add  remotely  sensed  SSH 
anomalies  and  subsurface  data.  The  goals 
of  the  NOAA  system  are  to  provide  (a) 
accurate  estimates  and  forecasts  of  the 
coastal  ocean,  (b)  initial  and  boundary 
conditions  to  NOAAs  regional  and  coast¬ 
al  models,  (c)  coupled  circulation-wave- 
storm-surge  models,  and  (d)  coupled 
atmosphere-ocean  hurricane  forecasts. 

The  generalized  coordinate  approach 
used  in  HYCOM  minimizes  the  liabilities 
associated  with  a  single  coordinate  sys¬ 
tem  and  provides  the  user  with  the  flex¬ 
ibility  to  tailor  the  model  to  the  specific 
application.  Generalized  hybrid  vertical 
coordinate  ocean  models  are  currendy 
used  for  an  increasingly  diverse  suite  of 
applications,  from  high -resolution  now- 
casting  and  short-term  prediction  of  the 
regional  ocean  state,  to  global  tidal  simu¬ 
lations,  El  Nino-Southern  Oscillation 
(ENSO)  forecasting,  multi-century  cli¬ 


mate  simulations,  and  theoretical  stud¬ 
ies  of  the  ocean's  dynamics.  U.S.  ocean 
modelers  have,  therefore,  engaged  in  a 
broader  dialogue  as  to  whether  the  vari¬ 
ous  community  modeling  efforts  could 
be  channeled  into  an  ocean  modeling 
environment  (see  white  paper  on  HOME 
[Hybrid  Ocean  Modeling  Environment] 
available  at  http://www.hycom.org/ 
publications.html).  A  modeling  environ¬ 
ment  is  defined  here  as  a  uniform  code 
comprising  a  diverse  collection  of  inter¬ 
changeable  algorithms  and  supporting 
software  from  which  particular  models 
design  can  be  selected  (e.g.,  HYCOM). 

It  would  not  only  provide  diversity  of 
modeling  approaches,  but  also  would 
standardize  coupling  with  other  models 
(e.g.,  atmospheres  or  ocean  sub-models) 
and  various  ways  for  user  specification  of 
parameters,  grids,  domains,  initial  condi¬ 
tions,  forcings,  and  diagnostics. 
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